Mapping of strongly correlated steady-state nonequilibrium to an effective equilibrium 
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By mapping steady-state nonequilibrium to an effective equilibrium, we formulate nonequilibrium 
problems within an equilibrium picture where we can apply existing equilibrium many-body tech- 
niques to steady-state electron transport problems. We study the analytic properties of many-body 
scattering states, reduce the boundary condition operator in a simple form and prove that this 
mapping is equivalent to the correct linear-response theory, fn an example of infinite- [/ Anderson 
impurity model, we approximately solve for the scattering state creation operators, based on which 
we derive the bias operator Y to construct the nonequilibrium ensemble in the form of the Boltz- 
mann factor e~'^^^~^' . The resulting Hamiltonian is solved by the non-crossing approximation. We 
obtain the Kondo anomaly conductance at zero bias, inelastic transport via the charge excitation 
on the quantum dot and significant inelastic current background over a wide range of bias. Finally, 
we propose a self-consistent algorithm of mapping general steady-state nonequilibrium. 

PACS numbers: 73.63.Kv, 72.10.Bg, 72.10.Di 



I. INTRODUCTION 

Quantum transport theorjiS*^ has recently been one 
of the most exciting fields of condensed matter physics. 
The field not only has a great promise for providing the- 
oretical framework for the design of nanoscale electronic 
devices, but also has been a frontier of quantum many- 
body theory in nonequihbrium. 

In regards to predicting behavior of nano-electronic 
devices the current theories based on non-equihbrium 
Green function technique^ii have been quite successful. 
However, as we approach strongly correlated regime of 
transport^'^. under extreme nonequilibrium conditions, 
approximate Green function techniques become increas- 
ingly unreliable. One can make an analogy to the equilib- 
rium quantum statistical mechanics, where the combina- 
tion of diagrammatic methods and computational tech- 
niques has been vital for the great success in the quantum 
many-body theory during the past decades. 

Therefore, it is a pressing issue that we have an alterna- 
tive theoretical scheme to implement the computational 
techniques to nonequilibrium problems. Recently new 
approaches have been proposed to solve nonequilibrium. 
Mehta and Andrei^ have applied the non-equilibrium 
Bethe Ansatz technique to solve the interacting resonant- 
level model. Renormalization group theory has been 
extended to nonequilibrium systems^ to investigate the 
time-dependence^ . 

One of the main problems in applying numerical tech- 
niques to transport problems has been the lack of under- 
standing of nonequilibrium ensemble, specifically, as to 
how nonequilibrium boundary conditions can be imple- 
mented as a statistical operator. In steady-state nonequi- 
librium, Zubarev^ has extended the Gibbsian statistical 
mechanics to incorporate steady-state boundary condi- 
tions in the density matrix formalism. Later, Hershfieldi^ 
has shown that the nonequilibrium ensemble can be ex- 
pressed in the effective Boltzmann factor e~^^^~'^'^ with 
the bias operator, Y , in terms of the scattering state op- 



erators. Once the effective Hamiltonian H — Y inside the 
Boltzmann factor is obtained, one can use the equilib- 
rium numerical techniques to sample the nonequilibrium 
ensemble. 

Despite the great potential, application of the method 
has been limited due to the difficulties in constructing 
the bias operator—. Recently, the method has been im- 
plemented by the author^ in the system of quantum dot 
(QD) coupled with electron-phonon interaction. In the 
work, he has approximately constructed the bias opera- 
tor Y by truncating the non-local interaction arising from 
the nonequilibrium boundary condition and by expand- 
ing the scattering state operators up to the harmonic or- 
der of the electron-phonon interaction. However, in the 
strongly correlated limit, such perturbative expansion of 
the scattering state may not be sufficient to reproduce 
the many-particle transport such as the Kondo anomaly 
conductance. 

Ultimately, it is desirable that we develop a general 
computational scheme of mapping nonequilibrium sys- 
tem to an effective equilibrium Hamiltonian, so that equi- 
librium techniques such as quantum Monte Carlo method 
can be used to directly calculate nonequilibrium proper- 
ties. However, to achieve the goal, we must have a good 
understanding of the properties of the nonequilibrium en- 
semble. Therefore, in order to gain physical insights of 
the mapping, we first study the diagrammatic construc- 
tion of the nonequilibrium ensemble. 

The first main goal of this paper is to derive general 
analytical properties of scattering state operators, espe- 
cially the anti-commutation and completeness relations 
in the interacting limit, and derive various expressions 
of the boundary condition operator. By clarifying the 
fundamental questions on the properties of scattering 
state operators, we lay a foundation of the mapping the 
nonequilibrium. Furthermore, we show that the nonequi- 
librium mapping produces the correct theory in the zero- 
bias limit. 

The second goal is, by taking an example, to address 



what correlation effects have to be included in the bias 
operator Y to properly describe the strongly correlated 
transport. To this end, we consider the Anderson impu- 
rity Hamiltonian as a model for Kondo dot systems. We 
might expect that, since the original Hamiltonian has a 
strong on-site interaction, the strong correlation effects 
will be described simply by shifting chemical potentials 
in the source, drain reservoirs. However, it will be shown 
in this work that it is essential to include the correla- 
tion effects at the level of Hamiltonian in the boundary 
condition to produce the Kondo anomaljii^. 

This paper is organized as follows. In Sectional we 
introduce the idea of mapping the nonequilibrium and 
present the analytic properties of the scattering state and 
the bias operators. The proofs are provided in Appen- 
dices. In the following Section HTll -A. we develop a pro- 
cedure to expand the scattering state operator using the 
slave-boson representation of the QD states in the infi- 
nite on-site Coulomb interaction. We derive expression 
for the scattering state operators and give physical inter- 
pretations in the expansion. In the next subsection llllt 
B, we construct the bias operator Y and explain how 
the nonequilibrium Hamiltonian is calculated in the non- 
crossing approximation in the equilibrium picture. In 
the subsection Illlt C. we derive an approximate expres- 
sion for the current. In section Hvl we present numerical 
results and discuss strong correlation effects in the dif- 
ferential conductance. In section we discuss how the 
limitations in the previous section can be improved in a 
self-consistent algorithm. 




(b) 



FIG. 1: (a) Schematic energy diagram with chemical poten- 
tial bias <E> applied between the source and drain reservoirs. 
The empty level of the quantum dot is aligned at the mid- 
dle of the two chemical potentials, (b) After the mapping of 
steady-state nonequilibrium to an effective equilibrium via the 
boundary condition operator Y , the energy diagram oi H — Y 
has equal chemical potential. The non-zero current is recov- 
ered in the modified hybridization functions (solid curves) be- 
tween the QD and the reservoirs. 



II. ANALYTIC PROPERTIES OF THE 
MAPPING 

In this section, we discuss the general properties of 
scattering state operators with many-body interactions 
present in quantum dot systems and obtain expressions 
for the bias operator in a simple form. We prove that 
the zero-bias limit in the mapping of nonequilibrium pro- 
duces the Kubo formula. The detailed derivations are 
presented in Appendix. 

The general scattering theory has been well-studied 
and its formalism can be found in many textbookai^. 
However, most of them are studied in terms of state vec- 
tor notation and here we present some of the relations 
in the operator form. We first consider a quite general 
model of quantum dot systems. The constraint is that 
the many-body interaction is confined to the quantum 
dot. The Lippmann-Schwinger scattering state operator 
^aka ^^ expressedi^ in the operator equation. 



•0^ J. ^ c} , + 
^ak<T aha 



1 



Eafc 



C + irj 



^Vclka^ 



(1) 



where 



-'aka 



is the free continuum state creation opera- 



tor inside the reservoir a at the asymptotic energy eak 
with the continuum index k and spin a. ry is a infinitesi- 
mal positive number to define the asymptotic limit of the 



scattering state. ^^ Here the Liouville operators for the 
full Hamiltonian C and the interaction Cv are defined in 
the commutation relation such as £A ~ [H, A] for any 
operator. The a{= ±) index refers to the reservoir [+ ^-s- 
source reservoir (L), — <-> drain reservoir (i?)]. 

Since the Hamiltonian has many-body interactions, 
the scattering state operator tp'' has terms with many- 
particle operators, and it is not clear if they obey the 
regular fermion anti-commutation relation. In non- 
interacting models, it can be shown straightforwardly 
that the anti-commutation relation holds. It is shown 
in Appendix A that the following commutation relation 
holds even in the interacting limit: 



iV'lfca' V'a'fcv} = Saa'Skk'Sa 



(2) 



Therefore, ^p',ip work like typical fermion operators and 
any scattering state created by repeated application of 
i/jt's on a particle- vacuum can be considered a many- 
body scattering state. 

One of the obstacles in implementing the mapping of 
nonequilibrium has been the lack of understanding of the 
bias condition operator suggested by Hershfield^°, 



Y 



$ 



2^ 



(*L 



Lka^Lka " ^Rka'^Rka 



(3) 



We reduce the Y operator in a physically appealing form 



in Appendix B as 



y = $ 



E- 



akcr 



^Oiko 



-c 



IT] 



(4) 



they contain many-body interactions to infinite orders in 
the original basis. However, after the summation over all 
(aka) indices, only finite order terms up to the original 
Hamiltonian survive. 



with the current operator / which will be defined later. 
Although the bias operator has been derived from the 
boundary condition imposed with respect to the chem- 
ical potential difference, the above formula relates the 
potential-driven ensemble to a current-driven ensem- 
blei^. 

Expressing the bias operator in terms of current oper- 
ator becomes particularly useful in proving the zero-bias 
limit of the conductance. The formulation of nonequi- 
librium mapping via Y must produce the same Kubo 
formula in equilibrium theory. The proof in Appendix C 
shows that the mapping of nonequilibrium has the correct 
description of transport physics in the low-bias limit. So 
far, the present formulation of nonequilibrium has been 
shown to be correct in the non-interacting models at fi- 
nite bias^^ and in the interacting models in the zero-bias 
limit. 

The notion that we can treat the scattering states as 
independent dynamical degrees of freedom has been the 
central idea in the mapping of nonequilibrium. In the 
non-interacting modeU^, the scattering state operators 
are on a firm ground due to the anti-commutation rela- 
tion mentioned above and the completeness of the scat- 
tering states encapsulated in the following identity 



]'iplkai'o.ka 



/ . '^aka'^aka 
aka 



Y^dUo, (5) 



provided that there exist no bound states, which is the 
case in the limit of large bandwidth. In the interact- 
ing systems, the scattering state operators lAlfctr become 
dressed with multi-particle scattering. Therefore, i/'Ifca 
should be defined in the many-particle basis and the va- 
lidity of the above equation is not obvious. In Appendix 
D, we derive the relation for a fairly general class of in- 
teractions. 

Similarly, behind the idea of writing the bias operator 
Y by shifting the chemical potential of ilr^j.^ , we expect 
for a consistent theory that 



H -Y 



E(^ 



^ka - a$/2) V'LcrV'afco 



and therefore 



-H" = ^ eakatpikcr'^aka 



(6) 



(7) 



aka 



We derive a general relation for the above summation 
in Appendix D and explicitly demonstrate that the re- 
lation Eq. (UI) holds for the Anderson impurity model. 
This relation corresponds to the intertwining relation in 
the S'-matrix formalism in the absence of bound states. 
Even if the products '(jy^f.^'ipaka has a quadratic form. 



III. NONEQUILIBRIUM IN ANDERSON 
IMPURITY MODEL 



For the remainder of this paper until Section we 
will consider an example of electron transport in the 
infinite-C/ Anderson model as a model for strongly cor- 
related transport system. In such systems where lo- 
cal interaction dominates the electronic structure, it is 
natural to choose the local states as the unperturbed 
states. One of the very successful and popular meth- 
ods is the slave-boson technique, often used in the limit 
of infinite interaction strength. In our model, by letting 
the on-site Coulomb interaction to be infinite, we effec- 
tively project out the Hilbert space associated with d? 
QD states. Therefore the available local Hilbert space is 
{|0), jdcr) : a = 1, • • • ,iV}, constrained by the relation 



N 



|o)(o| + ^M.)(d. 



= 1. 



(8) 



Due to the truncation of the Hilbert space, the QD states 
I do-) cannot be represented by ordinary fermion oper- 
ators. One of the techniques to overcome this prob- 
lem is the slave-boson method, where we associate the 
empty state |0) by a bosonic state h^\vac)^ created by a 
slave-boson creation operator w on an imaginary vacuum 
\vac) and the singly-occupied states \da) by fl\vac) with 
pseudo- fermion operator /J. The above constraint can 
be now recast in the following operator relation. 



N 



b'^b 



/ , J a J a 



(9) 



Under the assumption of the above projection, the Hamil- 
tonian can be written in terms of fermion-boson opera- 
tors as 

H = Ho + V (10) 

Ho = ^eafcc|,fe^Cafc„-|-ed^/]/^ (11) 



aka 



aka 



^ - E^(^/^^"^- + ^^.^/'^^' 



(12) 



with the QD level energy e^, tunneling amplitude taka- 
and the volume factor f2. The unperturbed Hamiltonian 
Hq has the reservoir continua and local basis as discon- 
nected. The perturbation term V turns on the tunneling 
between the reservoirs and the quantum dot. 

Since the steady-state ensemble must be time- 
independent, we construct the nonequilibrium boundary 



condition in terms of the scattering statesiS. The scat- 
tering state creation operator f/'lfco- '^^^ ^^ written in the 
Lippmann-Schwinger equatior&iSiii: 






1 



k(T 



aka 



^ak - Co^if] 



-^v€k.^ (13) 



with the LiouviUe operator Cq with respect to Hq de- 
fined as CqA = [Hq, A] for an arbitrary operator A, and 
similarly CyA = [V,A]. This formula is equivalent to 
Eq. ^. The scattering state operator is expanded in a 
perturbation series of the tunneling part Cy , 





Vaka — ^aka + /^ Vn,ak(T 
n=l 
oc y -. \ n 



The n-th order term can be derived recursively as 

1 



Vn,ak 



£ak - Co+iri 



^Vfpl-l,aka- 



(15) 



The statistical operator Y for the nonequilibrium bound- 
ary condition is defined as in Eq. ||2J). 

If the interaction is local, one can easily see that 
the scattered part of the wave- function, A-0^j,^ = 

J2n'^i,aka' ^^^ ^^^ fc-depcudence only through e^fc. 
To eliminate unnecessary band-edge effects and bound 
states, we consider a infinite band Lorentzian density 
of states (DOS)^* for the both reservoirs Na{e) — 
{D/iT){e^ + D2)-i with the half-bandwidth D. Then 
at fife — ^Rk' we have A-f/'j^^.j^ ~ A'0}jj,,jj. After the k- 
summation with the DOS of broad bandwidth and by as- 
suming the symmetric source-drain parameters, we have 
a cancellation between Aip'j^i^^AipLkcr and A^/'jjj.^A-f/'fl.fccr 
inF: 

^ = ^0 + ^J2 {^^Ika^^Lka " Mnk'aCRk'a + h.C. 
k(T 



with fo = ($/2) Ekai^^lkaCLka " C^RkaCRka)- 



(16) 



(d) 



(e) 



aka fa a'k'a }u aka h a'k'a' b 




akcj fa b fa 



ka b fa' 



FIG. 2; Diagrams in the expansion of the scattering state 
operator V'Ifco-' expanded perturbatively with respect to the 
tunneling term V in the infinite-[/ Anderson impurity model. 



and the tunneling amplitudes are independent of fc, i.e., 

tq=t. 

In the n = 2 term, we first compute Cv{bf^) as 

q' a' 

where b^ ~ Q has been used in the above calculation since 
any physical state subject to the constraint Eq. is pro- 
jected out by the application of b^ operator. Therefore, 



A. Expansion of tpl^^.^ 



/,t 



We expand the scattering state operator "^afco- as de- 
picted in FIG. 13 FIG. El^a) shows the diagram for the 
n = 1 contribution. The solid line is for the conduction 
electron, the dashed line for the pseudo-fermion f^ and 
the wavy line for the slave-boson b. 



Yl.aktT ~ 



1 



eak - Co + irj 



^vclt,^ 



^a. 



bfl 



Vncak - ed + iri' 
(17) 

From now on, we use a shorthand notation q = (ak). 
As mentioned above, we assume the left-right symmetry 



^iqa 



ti>ci,,,{s,,,+b^bs,,,-un) 

' +*'7)(eq -f-d + ivi) ' 






(19) 



The second term in the parenthesis of the numerator cor- 
responds to the diagram in FIG.[2Ib) and the third term 
to (c). 

The third order term {n = 3) is obtained from -03 = 

(eafe — £o + "?)~^'CyV'2,9cr- We have simpHfied the algebra 
of commutation relations regarding the QD operators by 
imposing the constraint Eq. lO, such as bb^b = b and 

bfl'fa'fl = bflScrcr'- With regards to the conduction 
electron operators, we ignore the off-diagonal density and 



electron (hole) pair operators as 



Cq"(T"Cg,^, 

t t 



Cq'a'cl,„,Sq'q"Sa>cr" (20) 



0. 



(21) 



If we do not consider any effects on transport from 
spontaneous magnetic ordering in the continuum states, 
^•^- ('^Qfea'^a'fe' -cr) ~ '-'' '^^ pairing correlation, i.e. 

\ a 



ka a'k'a' 



0, it is reasonable to drop the correspond- 
ing operator products. 

Justification of the q' — q" condition seems less ob- 
vious. We discuss this in conjunction with the U = oo 
conditioniS as follows. As will be discussed below, we 
will replace the operator product c ,^,Cq"a' by its expec- 
tation value, which amounts to contracting the out-going 
electron- hole lines in FIG. |21 (f-g). However, if q' ^ q" 
there must be a slave-boson-fermion loop to mediate the 
off-diagonal conduction electron states. This contributes 
to an extra-charge Q number in the Coleman's formu- 
lationiS of the infinite-C/ model and can be projected 
out. It can be argued that the outgoing conduction elec- 
tron line will be eventually contracted in higher order 
diagrams. However, that corresponds to the so-called 
crossing-diagramsSS which are of higher order in 1/N and 
negligible in the large N limit. Therefore, in the context 
of the large-[/ and large- A^ limit, the off-diagonal density 
operators are ignored. 

With above approximations, the third order contribu- 
tion of the scattering state operator becomes 



^L-^ = 



tc./vTi 



(fg ~ £d + iv^ 



-zr + S,(e,))6/t, (22) 



where T ^ Tl + Tr, T^ = Q-^Trt^J^k^i^ ~ ^c.k) = 
TTt'^N{0) with iV(0) the density of states of the reservoirs. 
The retarded self-energy operator is defined as 



n ,^ eq - tq' + IT] 



(23) 



Applying the same approximations discussed above, we 
obtain higher order perturbation terms in a geometric 
series as 



V'odd, 



E ^: 



t 

n,q(T 



(24) 



n—1,3.--- ,cx) 



^bfl(ec,k-ed + tr-±r{eq)) '.(25) 



go- 



So far the self-energy operator is derived from the equa- 
tion of motion and is independent of the boundary con- 
dition. 

Here we replace the pairs of electron operators in the 
numerator of the self-energy operator with average val- 



ues in the asymptotic Hmit, cj ^c, 



q,cT^q<^ 



(V'L.xV'afeo 



/(cafc — a$/2) = fa{eak), with the Fermi-Dirac func- 
tion /(e) = (1-1- e^'^)~^ . The different chemical potential 



of the source-drain reservoirs is taken into account in the 
a<I'/2 term. By taking the expectation value over the self- 
energy operator, we have incorporated the boundary con- 
dition in the scattering state. We define an effective re- 
tarded Green function gdi^) = [^ak — £d + i^ — S''(eq)] 
in terms of the retarded self-energy S'^(e) given as 



^'■(.g) = ^ E 



tlU'ieq') 



^ ,^ ''q 



ITj 



, 1^ J eq - e' + iri 



(26) 



with the renormalized DOS for the a-th reservoir 
ria'ie') = Na'{e')/Na'{0). The self-energy E''(e) has 
a logarithmic singularity in the limit of T = for 



a<i>/2 as 



de 



,n^'{e')fc'{e') 



In- 



IT] 



a$/2| 



iTTna'{eq)9 {a^/2 - 



(27) 



This has the familiar form of the Brillouine-Wigner per- 
turbation theorjiSii and goes back to the typical expres- 
sion in the zero-bias limit (with a single chemical poten- 
tial) . We now write the odd-order perturbation terms of 
the scattering state operator as 



A^l 



dd.qa 



9d{eq)bfl, 



(28) 



which has the same form as in the non-interacting reso- 
nant level QD system when gdi^q) is replaced with the 
non-interacting function^^ (cafe — Cd + *r)~"'^- The renor- 
malization factor due to S'^(e) is essential to construct 
the bias operator Y in the strongly correlated limit. As 
will be shown later, anomalous Kondo peak at zero bias 
cannot be produced without this logarithmic singularity. 

The approximations used here are motivated to pro- 
duce physically intuitive formula, similar to the scat- 
tering states in the non- interacting limit. We view 
that the many-body particle terms in the scattering 
state operators as out-going single particle states dressed 
with contractible many-particle excitations. One can in- 
terpret the replacement of many-particle operators by 
single-particle operators with renormalized amplitude as 
a quasi-particle approximation, in the similar spirit of the 
Fermi-liquid theory. 

The even order perturbation terms can be obtained 

from AV-lvon.ga = i<^q-'^0 + 'iv)'^^V'^ijlM,qa ^^^' ^siug 

the approximations leading to Eq. (|28|l . we get 



^2 

T^cvcn,g(T Q / J 



i(eg) 



Ql J 



Q ^ e„ — £„' -I- in 

q'lj' ^ ^ 



4,^,(5,,,+5t6J,,,-/,,/t) 

(29) 

which can be interpreted as the renormalized scattered 
wave. Finally the approximate scattering state operator 






^9d{e,)bfl 



n 



E 

q'a' 



9d{eq) 

- f-q' + i-q 



(30) 



which has the similar structure of the non-interacting 
scattering state. 

The physical meaning of the above scattering state is 
quite transparent. The second term is for the conduction 
electrons to tunnel on to the QD site with the amplitude 
given by the full Green function. The last term involving 
continuum states is cither from the potential scattering 
or from the exchange of the conduction with the QD 
electrons. Due to the approximation of [/ = cxd, local 
electrons first tunnel out to the reservoir and the empty 
local state is subsequently filled by a conduction electron 
from the reservoir {d} -^ dP —* d}). 

If the outgoing conduction electron were to be in the 
same spin state as the incoming spin, there are two pos- 
sibilities. First, when the initial QD is in the empty 
state b^\vac), a conduction electron goes through the 
potential scattering to hop onto the QD and then tun- 
nels out with the same spin (d" ^ d^ ^ d°). It can 
be seen explicitly from Eq. Ij^HI that the initially empty 
QD state (5^|vac)) does not change after the scattering 
[{5aa' + fe^Wcrcr' - fa'.fl)b^\vac) = 6aa'b'^\vac)]. Secoud, 
if the QD was singly occupied with the same spin as 
the incoming state (/J|uac)), the QD electron exchanges 
with the incoming electron without flipping the spin, 
[(1 + 6^5 - f,fl)fl\vac) = fl\vac)]. Therefore, when 
the incoming and out-going waves have the same spin, 
the factor [5cru' + b^bSaa' — fa' fa) can be considered as 
one. 

However, if the QD had a different spin state (/^, |wac), 
cr' 7^ a) from the incoming electron, the outgoing electron 
cannot be in the same spin state as the incoming one 
[{1 + b^b— fcrf^)fl,\vac) ~ 0], since ^^-configurations are 
forbidden by the U — oo condition and therefore the QD 
electron must hop out first in a scattering event. 



B. Effective Hamiltonian and Non-crossing 
Approximation 



Caka'^a'k'a', Created by the steady-state boundary con- 
dition. Finally an approximate effective nonequilibrium 
Hamiltonian H = H — Y becomes 



n = H-Y 



E 



^cxk 



a$ 



'^aka'^akn + ^d 



/ , faJc 



t v^ 

+ TT^E 



y/n 



aka 



a$ 



1 -Z-gdi^ak) ) bflcakcr + h.C. 



(32) 



Yr, 



Since the scattering state Eq. H31I) is an extended state, 
the product of iplg-i'qa in Y produces non-local interac- 
tion between all possible basis states and Hamiltonians 
with only local interaction become non-local after the 
mapping. Usually we need a truncation scheme to treat 
the non-localitjii^. In the above Hamiltonian, the non- 
local effects are implicitly included in the modified tun- 
neling and we do not consider additional non-local effects 
in this work. 

In this section, we calculate the the on-site QD Green 
function Qdi^)- Since the non-local terms in Kcc have 
indirect effects to Gd{^), we temporarily drop Ycc in the 
calculation of Gdi^^)- We use curly symbols, such as TC 
and Gd{'-^), to denote the quantities which are derived 
with H — Y as the time-evolution generator. 

To solve the effective Hamiltonian subject to the con- 
straint condition Eq. @, we use the Coleman's projec- 
tion method^^ and perform the nonperturbative partial 
summation of Feynman diagrams in the non-crossing ap- 
proximation (NCA)2£. The constraint condition is im- 
posed within the grand-canonical ensemble scheme with 
respect to the quantum dot charge operator Q defined as 
Q = 6^& + J2a fafcr- Only the subspace Qi (with Q = 1) 
is physically meaningful. The grand-canonical ensemble 
is introduced with a fictitious chemical potential —A as- 
sociated with Q in the partition function Zq{X) as 



ZG(A)^Tre-««+^'3) = ^Z,(Q)e 



-I3XQ 



(33) 



Q=0 



with the partition fmiction Zc{Q) in the subspace of the 
fixed charge Q = Q. The expectation value of an opera- 
tor A averaged over Z^Q ~ 1) is accomplished by taking 
the A — > oo limit as 



Here, we focus on how the reservoir-QD tunneling is 
modified by the bias. The modified tunneling results 
from the second term in Eq. (|31|l . and the bias opera- 
tor Eq. H16() is written as 



i^ = ^2^a Clka'^aka + -7=gd{^ak)bflcaka + h.C. 



a.k 
+Ycc, 



Vn- 



(31) 



where the last term Y^c accounts for the non-diagonal 
coupling between conduction electrons, proportional to 



(i) = hm m^. 

^— {Q)x 



m 



where the average {■■■)\ is taken over the grand- 
canonical ensemble before A is taken to the infinity and 
regular Feynman diagram technique can be applied to 
evaluate the average. 

Feynman diagrams are partially summed by treating 
the inverse of the QD degeneracy, 1/iV, as the expansion 
parameter. The non-crossing approximation exploits the 
fact that diagrams with crossing one-particle propagators 
contribute to the high order expansion in 0{\/N'^), and 
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FIG. 3; Non-crossing approximation diagrams for the self- 
energies of slave-boson (6) and pseudo-fermion (/o-)- 



therefore can be ignored in the limit of large QD degener- 
acy. This procedure can be summed up in the self-energy 
diagrams depicted in FIG. 13 The dashed line represents 
the full Green function for the pseudo-fermion, the wavy 
line for the slave-boson, and the solid line for the bare 
conduction electron. After A — j oo limit is taken, the 
self-energies Sf,(w), S/(w) for h, f^ can be expressed as 

E,H = iV^ I dJ^^Ue)Gf{oj + e'^) (35) 



a=±' 



with e^ = e — a^/2. Gbiuj), Gf{uj) are full one-particle 
Green functions for the slave-boson and pseudo-fermion, 
respectively, defined as 






(37) 
(38) 



The effective hybridization function Ta{e) is defined with 
the modified tunneling amplitude in Eq. I|32(l as 



fa(e)=^i' 



a$ 



5rf(e) 



iV„(e). 



(39) 



The closed set of Eqs. (I35I39|I is self-consistently solved. 
The shifted continuum energy levels and the effective 
hybridization function in Ti. is schematically sketched in 
FIG. ^b). The dashed line is the original hybridization 
function at zero bias with a flat DOS. As the bias <i> 
is turned on, the hybridization is modified by the term 
proportional to ^gd{c) in Eq. H39() . as drawn in solid line. 
The enhancement of the effective hybridization from the 
source reservoir Tl{uj) at lower energy can be understood 
as follows. The bias voltage creates a current-carrying 
nonequilibrium ensemble. However, after adjusting the 
chemical potentials of the two reservoirs to the same 
level, the electron current (L -^ R) should be restored 
via the effective tunneling. By enhancing the tunneling 
magnitude, hence the hybridization, to the source reser- 
voir and suppressing the drain-hybridization for the low 
energy states, we recover the electron flow from L to R. 
Similarly, the enhanced drain-hybridization at higher en- 
ergy promotes the hole-current from Rio L. In fact, it is 
not only the magnitude of the tunneling parameters but 
also the phase factor (coming from the retarded Green 
function) associated with the tunneling electrons which 



contribute to the current. This phase factor breaks the 
time-reversal symmetry. 

The on-site Green function of the QD state, Qd{Lo), can 
be expressed as a convoluted integral of 6 — / spectral 
functions as 



Gd{(^) = / de- 
J t 



Pdie) 



UJ ~ e + i?] 
with its spectral function pd(e) given as 
1 -I- e-'^' 



(40) 



Pdie) = 



Z 



f 



de'p,{e')pf{e + e'), (41) 



where Pb{(), Pf{^) are spectral functions of the 6, / oper- 
ators, respectively, and the local partition function Z f is 
given as 

Zf = Jdee-^'[pb{e) + Npf{e)]. (42) 



C. Cur rent- voltage relation 

The current operator can be derived from the conti- 
nuity equation d{ed'ld„) / dt + Ij^^ + Iji^ — 0. The cur- 
rent operator through the left- and right-reservoirs are 
obtained as 

^ k 

= ^T.i^lk.b^f-^^f^'^'o.ka). (43) 

^ k 

In the steady-state condition, {d]^d(j) is constant and 

Il+Ir = (44) 

with la = X]<t(^"o")- I^ ^^ ignore the Ycc term in the ef- 
fective Hamiltonian as in the previous section, we violate 
the symmetry condition Eq. 144(1 as follows. Without the 
Ycc term in TL, the expectation value of operators appear- 
ing in Eq. (|43|l can be expressed in the finite-temperature 
Green functions as 



^ k uin k 



Gd{iL^n 



e„fe-Ha$/2' 

(45) 

where the QD-Green function Eq. (|40|l has been analyti- 
cally continued to the imaginary Matsubara frequencies, 
iujn = i{'2n + l)7r//3 at integer n. The on-site Green func- 
tion Gdiii^n) is independent of the reservoir index a. As 
discussed in the previous section, the tunneling from the 
source reservoir is dominated by the electron-hopping, 
d^ — > d^, while the tunneling from the drain reservoir 
is dominated by the hole- hopping, d^ — *■ d°. There- 
fore, the modified tunneling in the source and drain- 
hybridization do not have the same values of at a given 
energy (See FIG.QJ and the expectation value computed 
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FIG. 4: Spectral function of the retarded Green function for 
the quantum dot site and the hybridization function from the 
left- and right-reservoir. The bias $ is set to 1, ti, = tn = 0.48 
and /3 = 1/T = 160. The Fermi energy positions for the left 
(w = "l?/2) and the right {uj = — $/2) reservoirs are marked 
by dashed lines. 



from Eq. H45|l does not satisfy the charge conservation, 
Eq. 1331). 

Physically, it is no surprise that the Ycc term is crucial 
for the current. Ycc has the direct coupling between the 
source and drain reservoirs driven by the bias and has the 
effect of closing the circuit. Furthermore, phase factors 
in the coupling act as current battery which forces a net 
current from the source to drain through the quantum 
dot. 

To restore the L ~ R symmetry, off-diagonal Green 
functions for conduction electrons G^ , (ia;„), resulting 

from Ycc, has to be considered in Eq. (|45|1 . 
it 



^^X!^^-^'^"^"'^'^^ 



a/0 



= tt^ 



rJ2oY1 ^"i' ii^n)Gdiii^n) 



(46) 



qq' 



Here the conduction electron Green function G^ ,{iujn) 
is a propagator in the so-called 'cavity' HamiltonianSi, 
where the QD site and its tunneling to the reservoirs 
are removed from the full Hamiltonian Ti. It should be 
emphasized that since the bias operator Y is built upon 
the correlated scattering states, even the 'cavity' part of 
the Hamiltonian has correlation effects. 

Recalling that much of the correlations effects are al- 
ready present in the full Green function gd(e) in the coef- 
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FIG. 5: Coupling between the continuum states in i? — Y" is 
reconstructed by introducing a QD cluster and fitting the QD 
Green function Qdij^)- 



ficient of the scattered wave [third term in Eq. %M\ \ , we 
approximate Eq (|31|l . in analogy to the non- interacting 
scattering state, as 



%. 



^qa 



q 



tata'gdi^q) f 



Eq' + iri 



(47) 



which amount to replacing the factor {Sa-cr' + h'^hSaa' — 
fa' fa) by Saa'- As discusscd at the end of Section ITTll A. 
{Saa' + b'^bSaa' — fa'fl) = Saa' without Spin- flip Scatter- 
ing. 

We now obtain a simple expression for the off-diagonal 
Green function G'^ ,{iujn) over the 'cavity' Hamiltonian 
based on the non-interacting scattering state operators 
V'Jcr using the following trick. The Green function gdi^) is 
fitted by introducing an arbitrary, but finitc^^ number of 
discrete levels which couple to the QD site [see FIG.|3Ja)]. 
For example, we prepare a finite cluster QD of arbitrary 
geometry, with the resulting energy poles distributed un- 
der the envelope function given by the spectral function 
of gd{Lu), [= —TT^Hmgdiuj)], as depicted in FIG.E^b). We 
choose the finite cluster such that the peak positions and 
their weight match the quasi-particle energies and the 
renormalization factor. Broad peaks can be fitted by in- 
troducing multiple cluster peaks underneath, as sketched 
in FIG.ISJb). The following derivation does not depend 
on a particular fitting scheme as long as the fitting can 
be made accurate. 

Then, the effective Hamiltonian 



h-y=y: 



^aka 






^lai') 



qa T^qcr i 



(48) 



can approximate the 'cavity' part Ycc in the original 
Hamiltonian Eq. ^Q. Although the QD parts of H -Y 
and H — Y can be quite different, it does not matter as 
long as we are only concerned with the 'cavity' Green 
functions. Carrying out the similar calculations as de- 
tailed in the Appendix of Ref»i^, the c—d Green function 
defined as 



GLrf(iw„) 



(49) 



= {CL. 



Ch- 



-S 



dt- 



C 



-CLc 



H-Y 



with clc = Vrj ^ Y.k CLk<y, is given by 
tL ^ 

[gd{ii^n + $/2) - gd{ii^n ■ 



= 2i 



GLd{i^n) ~ GdL{i^n) 



r 



$, 



The symmetry between the left and right is no'v 
Similarly, the QD Green function can be expre 

GdiiuJn) = 2 [5d(«w« + */2) + 9d{i(^,i - */2 

The 'cavity' Green function then satisfies tl 
equation 

tL GLd{i^n) - Gdhii^n) 

= tbtR [G\ji{iujn) - G?j.x(iw„)] GdiiuJn) 

TlTb. 9d{ioJn + $/2) - 9d{i(^n - <^/2) ~ 
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= 4i- 



r gd{ioJn + $/2) + .gd(iw„ - $/2) 



Gd(icj„) 



Finally, after replacing Gdiii^n) in the above equation by 
the nonequilibrium Green function Gd(i^n), analytically 
continued from Eq. H40|l . the Dyson equation for the full 
off-diagonal Green functions GLd{i^n) and GdL{i^n) be- 
comes 



= 4i 



tL [GLd{iUJn) ~ GdLii^n)] 

'^L^R 9diiuJn + */2) - .gd(icj„ 



(52) 



$/2) 



Qd{i^n 



From (c^^do-) 
expressed as 



r gdiiojn + «'/2) -f 9d{i^n - $/2) " 

~/?~^ X^n Gdhii^n), the current can be 



/($) 






4JVe FLFfl 

/i r/3 



E 



gd{i^n 



$/2) - gd{iuJn 



$/2) 



5d(«w„ + $/2) -I- .gd(it^„ - $/2) 



(53) 

^d(Jw„). 



The numerator gd{iuJn + *&/2) — gd{ii-^n — *&/2) acts as an 
energy window opened between the source-drain chemical 
potential ±$/2 and the Matsubara summation performs 
the integration over the window. This expression be- 
comes exact^^ in the non-interacting limit with the non- 
interacting retarded Green function g'^{z) as 

^W = ^^^E ^[■9S(*^«+*/2)-.9S(*^n-<i>/2)]. 



F P 



(54) 



IV. RESULTS 



We first discuss the retarded Green function of the 
quantum dot, (7d(e), which appear in the expression for 
the scattering state operator i^l^f.^, Eq. (J3IJ. FIG.^a) 
shows the spectral function, —Tr~^liiigd{LL'), at the unit 
bias <I> = 1 with tunneling parameters t^ = tn = 0.48 



FIG. 6: Spectral function of the thermal Green function with 
H — Y as the time-evolution generator. The bias $ is set to 
1, tL^tR = 0.48 and /3 = l/T = 160. 



and the inverse temperature (3 = l/T = 1/160. Through- 
out the paper, the half-bandwidth D of the Lorentzian 
DOS of the reservoirs is fixed at D — 4. Singular peaks, 
developed at frequencies w = $/2,— $/2, correspond to 
the Fermi levels of the left and right reservoirs, respec- 
tively. As the chemical potential of the drain is lowered 
by $/2, the ionization energy of removing one electron 
from the QD {E = Cd) to the Fermi level of the drain 
{E = -$/2) is reduced to A^ = \ed\ - $/2. Therefore 
a strong Kondo resonance forms on the electrons com- 
ing from the right-hand-side continuum at w = — $/2. 
Similarly, the resonance formed from the left reservoir 
produces a much suppressed resonance at lu = <i>/2. 

The resulting hybridization Eq. 1)39(1 due to the left 
and right reservoirs is plotted in FIG.E^b). As discussed 
below Eq. H39() , the hybridization to the source {L) reser- 
voir has strong intensity at low energy with a strong peak 
at w = — $/2 while it is suppressed at high energy with 
depleted intensity at uj = -|-<i>/2. The hybridization to 
the drain reservoir shows the opposite behavior. 

The I — V curves for the parameter set A^ = 5, e^ = 
"I) t^jj — 0.48 in the Kondo regime are shown in 
FIG. [TJa) at different inverse temperatures l/T. For 
all the tested parameter sets, the current monotonically 
increases with the bias. The differential conductance 
dl/d^ in (b) shows three distinct features: Kondo peak, 
inelastic current peak at <i> = \ed\ and a very broad in- 
elastic current at high bias. 

The Kondo peak, often referred to as the anoma- 
lous conductance peak, is the hall-mark of strongly 
correlated transport. We emphasize that the Kondo 
peak arises from the combination of the nonequilibrium 
Green function Qd(i^) and the retarded Green function 
gdiiu! ± $/2) with the strong correlation effects present 
in both Green functions [see Eq. H53|l ]. The conduc- 
tance G{^)[= edl/d^] is normalized to the unitary limit 
Go = Ne'^/h with the local degeneracy A^. 

To discuss the main conclusion and the limitations of 
this work, we present in FIG.|Hlcomparisons of the results 
with other approaches. First, to show that the strong 
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limit, the zero-bias conductance becomes 



G($ = 0) = 






(55) 



1 1.5 2 

Bias (4)) 



FIG. 7: (a) I — V curves (current in arbitrary unit) in the 
Kondo regime, (b) Conductance G{^) = edi /d^ normalized 
to the unitary limit Ne'^/h (with A'' the degeneracy in the 
Anderson model). Three peaks correspond to the Kondo res- 
onance, inelastic transport via the charge excitation of the QD 
at "1> = |ed|, and a broad background of inelastic transport of 
incoherent multiple charge excitations on the QD. 



correlation effects should be included at the level of the 
effective Hamiltonian H — Y, the logarithmic singular- 
ity in Eq. (|27|l is turned off by replacing the Fermi-Dirac 
factor fa{e) by one. The resulting I — V characteristics 
(dashed line in (a)) does not have the Kondo peak. It 
is clear that even though the Kondo peak is present in 
the nonequilibrium Green function Gdi^), its absence in 
the retarded Green function 5d(e) and therefore in Y, 
leads to no zero-bias conductance peak. This proves that 
the strong correlation effects in the bias operator Y are 
essential to build correct steady-state nonequilibrium en- 
semble. Furthermore, at large bias $ ~ 2|e/|, not only 
is the broad spectrum of inelastic transport missing but 
also the differential conductance becomes negative. This 
underestimation of current is due to that the ionization 
spectrum of the source reservoir at energy $/2 -I- e^ and 
that of the drain at energy — $/2 + e^ have reduced over- 
lap. However, proper inclusion of the retarded self-energy 
term compensates the reduced overlap through a strong 
amplitude in the ^^-configuration which mediates the in- 
elastic tunneling, (see the strong spectra at cj > for 
$ = 2|ed| in FIG.©. 

Since the Friedel-Langreth Fermi-liquid relation of the 
Anderson model has the QD spectral function at zero 
frequency p(0) = (I/ttF) sin^(7rnrf/A^)^^ in the infinite-t/ 



Using the typical value of n^ = 0.9 in these calcula- 
tions, we get G($ = O) w 0.29{Ne^/h), about 34% larger 
than the values shown in FIG. Even though this work 
shows that the correlation effects included in Y produces 
strongly correlated nonequilibrium ensemble, the discrep- 
ancy of the zero-bias conductance is problematic. This 
is due to the different approximations used to calculate 
separately the retarded Green function gd{^) in the scat- 
tering state operator Eq. (|26|l . and the nonequilibrium 
Green function Gd{'^) Eqs. H35I36|I within the NCA. Com- 
paring the sharpness of the Kondo peaks in the FIGs.01 
and El we can see that the correlation effects in the re- 
tarded self-energy Eq. I|26() have been significantly under- 
estimated. In evaluating the current, Eq. H53|l . the peaks 
in the two Green functions get convoluted in the Matsub- 
ara sum and the over-broadened peak in gdi^) has the 
effect of smearing out the sharp Kondo peak in Gd{^), 
leading to the underestimated zero-bias conductance. 

It should be emphasized that this shortcoming is due 
to the inconsistent approximations in the two Green func- 
tions, .grf(w) and tj^(a;), not the mapping of nonequilib- 
rium itself. If the same approximation had been used 
for the both Green functions, for instance by using the 
NCA summation scheme (see FIG.|31) for the self-energy 
in the scattering state operator (see FIG.|21), the zero-bias 
conductance should have been correct. As proved in Ap- 
pendix C, the mapping of nonequilibrium produces the 
correct limit of the zero-bias conductance. Furthermore, 
as will be discussed in Section we develop an algo- 
rithm with one Green function and the above problem 
will disappear. 

The central peak at $ = le^l in FIG.|7fb) is due to the 
inelastic transport via charge excitation. Analogous to 
the onc-phonon inelastic process in the electron-phonon 
coupled QDs"'^^, electrons tunneling from the source to 
the drain electrode lose energy of the voltage drop (~ 
$) while creating charge excitation on the QD with the 
ionization energy ~ |ed|. 

The spectral function (FIG. 01 at $ = has a Kondo 
peak near the chemical potential and the ionization en- 
ergy peak at iu — ed- If we naively expect that the 
current is proportional to the integral of the zero-bias 
spectral function between the source-drain chemical po- 
tentials, the single-charge ionization process would re- 
sult in an incorrect conductance peak at $ = 2|ec;|, 
since we model that the source/drain Fermi levels are 
displaced by ±$/2 with respect to the QD levels. The 
current calculated by assuming a rigid spectral function 
/,igid (X /dwpd(w,$ = 0)[/(w-$/2)-/(co + $/2)], plot- 
ted as thin lines in FIG.|HIb), does not show the inelastic 
peak at $ = le^l. 

The broad background of inelastic transport is due to 
the strong amplitude of dP configuration as noted above. 
At such high bias, the system is out of the Kondo regime, 
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FIG. 9: Normalized conductance to the unitary limit with and 
without the strong correlation effects in the bias operator Y . 
To turn off the correlation effect in Y, the logarithmic singu- 
larity is suppressed by replacing the Fermi-Dirac function in 
Eq. 12711 . The I — V curve without the correlation effects does 
not produce the Kondo peak in the conductance. 



as can be inferred from the strong spectral weight in the 
positive frequency {d^ -^ tf) at <& = 2|ed| in FIG.IHl De- 
spite the noise in the conductance due to the numerical 
differentiation at high bias, the trend of increasing inten- 
sity of the inelastic transport with the tunneling param- 
eter (see FIG. 13 agrees with the above observation. We 



note that the infinite- 1/ approximation becomes increas- 
ingly unreliable at high bias since the strong electron- 
channel of the current-influx to the QD {d^ — > d^) is 
precluded. 

Finally, FIG.IHlpresents differential conductance as the 
tunneling parameter is changed. As the tunneling in- 
creases, the Kondo resonance becomes stronger and re- 
sults in pronounced zero-bias anomaly peaks. Since the 
Kondo peak is shifted to the positive frequency from the 
Fermi energy by the Kondo temperature, the conduc- 
tance peak moves from zero to finite bias. 



V. SELF-CONSISTENT ALGORITHM 

It became clear from the above analysis that the strong 
correlation effects should be present in the boundary con- 
dition. The main problem in the above approach arose 
from the fact that the nonequilibrium ensemble and the 
time-evolution are governed by two different operators, 
H—Y and H. Therefore, to generalize the above method, 
it is essential to unify the two Hamiltonians with a con- 
sistent approximation to solve the generalized impurity 
model. In this section we propose a practical algorithm 
to solve this problem. 

We start from the observation that the scattering state 
operators, being eigen-operators of Hamiltonian, also 
serve as the basis for building the bias operator. The 
Lippmann-Schwinger equation Eq. (^ is equivalent to 
the equation of motion, 

(eafc-/: + i7?)V'L<T = «^c^fe<T- (56) 

From the anti-commutation relation Eq. Q, we have 



'CyV'L. = [tV'Lj = ^V'L 



(57) 



Combining the above two equations, we have a new equa- 
tion of motion 



i^ak 



a(f> 



- Ch-y + iv)^lka = i^clko 



(58) 



which is equivalent to the scattering state equation 
+ 1 



V'La 



^aka 



iak - a$/2 - Ch-y + iri 



r' T 

'-vC^aka- 



(59) 



We now have only one Green function, propagating with 
the effective Hamiltonian H — Y. The complication is the 
new tunneling interaction in Cy which is modified by the 
correlation terms in Y. However this effect becomes the 
second order effect since the energy denominator in the 
Green function has the main strong correlation effects 
in the modified hybridization function in Y, Eq. H39|l as 
considered in the previous section. Furthermore, in the 
small bias limit where the strong correlation effects are 
strongest, the correction in Cy is proportional to the bias 
$ and Cy « Cv becomes a good approximation. 

Therefore, we propose the following self-consistent al- 
gorithm. 
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(i) Start from an approximate Y, for example from the 
non-interacting resonant level model^^. 

(ii) Solve H — Y for the retarded QD Green function 

(iii) Approximate C'y — Cy and construct the new bias 
operator ip'^ using the quasi-particle approximation, 
Eq. (|31|) , and consequently Y with the modified hy- 
bridization function given by G'^f*'{Lu) [see Eq. (jSHJ]- 

(iv) Go back to the step (i) until it converges. 

In this scheme we have a single Green function and the 
problems arising from two Green functions, such as the 
underestimated zero-bias conductance, will not occur. 



through a chain immersed in a dielectric medium or 
conduction through fluctuating (classical) spin distribu- 
tions would be such examples. Since non-interacting 
scattering states can be written down exactly and the 
BoUzmann factors are now given in a definite expression 
exp[— /?(i/ — y)], we can study transport physics in such 
regime with the current level of formulation. 
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VI. CONCLUSION 



APPENDIX A: PROOF OF {ipt, ips} = <5r 



We have developed a procedure of mapping a strongly 
correlated steady-state nonequilibrium to an effective 
equilibrium. The anti-commutation and completeness re- 
lations of the scattering state operators are derived in the 
interacting limit and the bias operator has been rewritten 
in a physically appealing form. The mapping of nonequi- 
librium in the zero-bias limit has been shown to produce 
the correct linear-response theory. 

In the example of the infinite- C/ Anderson impurity 
Hamiltonian, as a model for Kondo dot systems, we have 
derived and analyzed the scattering state operators as 
the basis for the statistical operator which accounts for 
the nonequilibrium boundary condition. The strong cor- 
relation effects should be included at the Hamiltonian 
level in the statistical bias operator Y to produce the 
zero-bias Kondo anomaly. The current-voltage relation 
has been calculated using the equilibrium diagrammatic 
technique. 

Despite the growing understanding of the nonequilib- 
rium ensemble, there are issues to be resolved for this 
method to be readily applicable to general Hamiltoni- 
ans. To this end, we have developed a self-consistent 
algorithm with one kind of Green functions by exploiting 
the property of the scattering state operators as eigen- 
operators simultaneously to the Hamiltonian and the bias 
operators. 

The formulation of nonequilibrium steady-state 
through a mapping to an effective equilibrium not 
only provides an alternative way of understanding the 
nonequilibrium, but also can solve various boundary 
condition problems which are otherwise hard to formu- 
late. For instance, thermal boundary condition, such as 
thermo-electric effects under a finite temperature differ- 
ence between reservoirs, can be described quite naturally 
in the scattering state basis. 

One other interesting area is the electron transport 
with interactions with classical variables which are then 
self-consistently controlled by the transport, in the sense 
of the Felicov-Kimball modelHl. Electron transport 



In the following discussion, we will consider the Hamil- 
tonian with the interaction present only in the QD, such 
that Cycl = {t/^/n)d^ where Cy includes the tunneling 
and many-body interaction terms. For simplicity, we use 
the subscript s to denote collectively the reservoir and 
continuum indices (a, k) and suppress the spin index in 
the following. The Lippmann-Schwinger Eq. (Q becomes 
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t. 



1 



dt 



^/U er - C + irj 

ts 1 , 

^/Q es + C - iri 



(Al) 
(A2) 



where we have used the relation {CAy = {HA — AHy = 
A^H - HA^ = -CA^f. We write 



with 

rj l^Er ~ C + irj Eg + C — irj 
From the relations such as 

"1 



(A3) 




(A4) 



(cr - Eg - £ + iri) 



= d^Cs + ^ 



-(V 



1 



C + irj 
-rfM d. 



/Q \er — C + ir] 
Arfi can be written as 

1 f 1 



(A5) 



n 



d\d 



er — eg — C + iij \er — C + ii^ 
1 f .+ 1 



tr — tg — C + irj 
1 



d\ 



C + irj 



d\^ 



tg + C — itj 
1 



C — irj 



(A6) 
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Here we have used the relations {d^,Cs} — and 
£{Ab) = {£A)B + A{£B). The last term in the pre- 
vious equation can be rewritten as 



1 



-C 



-dl 



1 



-I 



€s + C- irj 



(A7) 



1 



er — C + 17] 



2ir] 



C — 17] 



which can be easily verified by multiplying (e^ — e^ — £ + 
2iri) to the both sides of the above equation. The first 
two terms in Eq. (|A6p are canceled by the third term in 
the 77 — > limit, A^s = Om^ Finally, we have the desired 
relation 



{lpl,1ps} = Sr 



(A8) 



which corresponds to the isometry of the M0ller operator 
in the S'-matrix formalismi^. 



Following the same steps as in the previous Appendix, 
we have 



V'LV'afc = cl^Cak + 



tn 



1 



VU-C + ir] 



(d^c^k-cl^d). (B5) 



Therefore the bias operator can be expressed in a simple 
form as 



y = $ 



yo 



1/2 



-C + IT] 



{d*CL — c'l^d — d'CR + cLd) 



(B6) 

with j/o == 5E/c(cLcLfc - cjjj^cflfc) and ^l = 

^^^^^^k^akCak- Defining the source-to-drain electric 
current as 



1 



{II - Ir) 



2 

= '-^{d^£L-cld-d^dR-cld), (B7) 

we have the bias operator in terms of the current operator 
as, 



APPENDIX B: BOUNDARY CONDITION Y IN 
TERMS OF CURRENT OPERATOR 

Following similar steps considered in deriving the com- 
mutation relation, we can reduce the bias operator Y into 
a convenient form in terms of current operators. We write 
y as 



$ 



Y = -{N'l- AAfl.), with Mc. = J2 ^ik^o^k- 



(Bl) 



Similar steps leading to Eq. I)A4|I give 



V'LV'afe = cli^Cak 






1 



tt 



1 



^ V ^ak — C + IT] 






tak + C-IT] 



(B2) 



The second and third terms can be rewritten as Eq. (|A5|) 
with er = €s = eak- 



1 



£ak- C + irj 
1 



-rft 



^ak 



--r--Utc„fc+ '°^7". dt-d ■ (B3) 
-£ + iJ] y e^k- L- + 17] J 

The second term in Eq. (|B3|I is canceled by the last term 
of the Eq. ||B2ll in the 77 ^ limit. 



1 

C-2ir] 



-dt 



17] 



1 



Cafe + 'C - 17] 



y = $j/ = $ 



1 



e —L 



IT] 



(B8) 



Or, alternatively, by making use of £v'(c]jj,CQfe) — 

(tak / V^){d' Cak ~Cakd) for Hamiltouians with only local 
interactions on the quantum dots. 



r = $ 



2/0 



-c 



— -Cvyo 
tr] 



(B9) 



which is equivalent to the equation of motion for the bias 
operator Y, [H,Y] — irj{Y — Yq), in Hershficld's paper—. 



APPENDIX C: DERIVATION OF ZERO-BIAS 
LIMIT CONDUCTANCE 

In the zero-bias limit, the present formulation via map- 
ping of nonequilibrium should recover the linear response 
formula with the conductance given by the current- 
current correlation function. With the expression derived 
in the previous section, one can prove this as follows. The 
current is evaluated via the thermal average taken over 
the nonequilibrium ensemble, 

/ = Z-'Tt (e-^(^-^)/) , with Z = Tre-P^"-^l 

(CI) 
By differentiating the current by the bias $ in y = ^y 
about the zero-bias limit, we obtain the expression of the 
conductance in terms of the imaginary-time correlation 
function, 



(B4) 



1 



eak - C + 17] 



dnd-d'f 



1 

eak + C - irj' 



-d 



81 f^ 

Go = e— - ey dT{y{T)I{Q)) with y(r) = e'^y, 



(C2) 
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where we have used the relation (/) = 0. 

The y term has two contributions, as in Eq. (|_B8|I . from 
yo and /. It can be shown that the contribution from yo 
can be dropped, as follows. For the sake of argument, 
we express the Hamiltonian H in terms of real numbers. 
Since the expectation values such as {yviScL) are all real 
numbers, they have no contributions in {yal)a according 
to Eq. IJB7|I . Therefore we have 



Gn 



dT(/l(T)/(0)), With/l 



-C + iri 



I. (C3) 



Analytically continuing the above integral from the 
imaginary-time to the real-time integral, we have the 
fluctuation-dissipation theorem 



dtx{t) with x{t) = -ie{-t){[i{Q),h{t)] 



Go = 

(C4) 
Fourier transform of the response function x(t) becomes 



x(w) 



dte''*+'"*/[/,e*^*/i]\ 



LO — C + i-q — £ -I- irj 



(C5) 



From 



LO — C + irj — £ + irj lu 
and 
Red 



-C + irj 



/ ) = Re 



— £ + 17] LO — C + irj 

(C6) 

i]i). (C7) 



-£ + irj 



the response function is written by the current-current 
response function Xii{t) = ^*^(0([-^(^))-^(0)]) as 



1 



Xiu;) = —xii{uj) 



1 

iuj 



1 



-£ 



-/ 



irj 



Finally we obtain the zero-bias conductance in terms of 
the linear response theory in agreement with the Kubo 
formula, 



lim 



X//(^) 



Go = limx(w) 



(C9) 



APPENDIX D: COMPLETENESS OF ^\ 

For non-interacting models^^ without bound states, 
the scattering state operators satisfy a completeness re- 
lation 



and the Hamiltonian is written as a sum over ipaka as 

X! ^"fcV'lfcaV'afeff = H. (D2) 

ak(T 

For the interacting case, it is not clear at all whether such 
relations hold when i/'Ifco- ^^ made not only of one-particle 
creation operators but also of many-particle excitations. 
First we assume that the Hamiltonian has many-body 
interaction H' term only on the QD (drr) site and it takes 
the form 



edY^dUo + H'{d\,d^,d\,di 



(D3) 



We derive corresponding relations for the interacting case 
starting from Eq. (|B5|) . After summing over the contin- 
uum variables on both sides of the equation. 



o/.k ak ak 

With the Hamiltonian Eq. (|D3|I . 

t, 



-£ + irj 



dUck ~ c^ferf<T 



(D4) 



C{dUa) = E ^ i^ikd. - dUck) + C'{dUa). (D5) 



ak 



Vn V 



Substituting this to the previous equation, we have 



E^'lfcV'afc = E^ifc"="fe 



1 



ak 



ak 



-C + irj 



= E^afcCafe+4c?a 



{-Cria- + C'Ua) 



£'n^(D6) 



ak 



-£ + irj 



= _ / dte''*+*"*(-^)/[/(0),/(t)]\. (C8) l^o) 
iuj ,1^, \ I the 1 



In the last line we have used that —C/{—C + irj) = 
l — irj/(—C + irj) —f 1 and assumed that there exist no iso- 
lated energy eigenstatesiSii^. If there is an isolated state 
Eq) outside the continuum it cannot be constructed by 
the scattering states. If ricr = d^d has a finite amplitude 
of ^5^0, with the creation operator rpl of the state \Eo), 
it does not contribute in the fc-summation of scattering 
states, i.e. —C/{—C + irj){rplrpQ) = 0. In other words, 
the operator irj/{~C + irj) projects out the contribution 
of isolated states in n^ . 

If the interaction is written in terms of d-density op- 
erators only, such as the Hubbard Coulomb interaction 
H' = UniTii, Clrir, = 0. Therefore, we recover Eq. HD1|) 
in the interacting case. 

We proceed similarly for Eq. 102^ : 



E ^afcV'LcrV'afc.T (D7) 

takl^ 



ak 



ak ak akcr ak 



15 



Since Cc^^^ = <^akc\ka + (iafe/V^)4> 
dl{CCaka) + {Cc^aka)'^<y = -<^ak{dlcaka ~ cl^.^da) , (D8) 

and 

dUCc^k.) + {Ccl,^)d, (D9) 

Therefore, Eq. (|D7p can be summarized, so far, as 



ykcr 



akcr 



/ . ^o^kC^f^^Cakcr + ^ /— (^"^(^q/cct + ^aka^' 



Oika 



ak<7 



+rzT7;;E^ [(^4Kfe. + 4,(/:rf.)] .(dio) 

akcr ^ 

From £4 = Eafe(*"'=/Vf^)Cafc. + ^^4 + -^'4, the last 
term in the previous equation becomes, before the spin 
summation is taken, 

_£^, Y. -^ [{Cdiy^ka + cikA^d,) 

ak 

[{Cdi)i-Cd^~edd^ + C'd„) (DU) 

4)iCd,)] 
[iCdi)iC'd,)-iC'di){Cd.)]. 



—C + irj 

+ {Cdl-eddl-C'di)iCd. 



By defining the non-interaction part as Hq we have 



^CafeV'L^V'afe-T 



(D12) 



^0 + TT^E [(^4)(>C'rf.) - {C'dt){£d.)] 



If the interaction is the Hubbard on-site Coulomb inter- 
action H' = Un^ng-, 

-{£dl){C'd^) + {£'dl){£d„) 
= U{Cdl)nsd„ + Unsdl{Cd,) 
= U [{CdDn^d, + 4 {Cin^d,) - iCng)d,}] 
= U [Cin^n^) - {Cng)n,] . (D13) 

Summing over spins, the interaction term becomes 

UY [(CdDiC'd^) - (C'dDiCd,)] = C{Un^ng) = CH'. 

(D14) 
Finally, we obtain for the Anderson impurity model 



E ''ak'lplkcr'^aka = Hf) + H' = H , 



(D15) 



-C + irj 



aka 



in the absence of bound states. 
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